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Abstract. We report the results of a broadband (0. 1-100 keV) 
X-ray observation of the nearby elliptical galaxy NGC 1052, 
performed with the BeppoSAX observatory. We confirm the 
presence of a bright (2-10 keV luminosity ^4 x 10 42 erg s _1 ) 
and strongly absorbed (/V H ~ 2 x 10 23 cm~ 2 ) X-ray source. 
The flatness of the X-ray spectrum (photon index, T, ~1.4), 
the estimated low accretion rate (m = M/M^dd ~ 10~ 4 ) and 
the radio-to-X-ray spectral energy distribution suggest that this 
observation may represent the first direct measurement above 
10 keV of an accretion-dominated flow in an elliptical galaxy. 
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1. Introduction 

The discovery of supermassive black holes (1O 6 -1O 9 M0) in 
the nuclei of several nearby galaxies (Magorrian et al. 1998) 
has raised the question of why most of them are not active. 
A possible solution is provided by the so called Advection 
Dominated Accretion Flows (ADAF; Rees 1982; Narayan & 
Yi 1995; Fabian & Rees 1995). In this class of accretion so- 
lutions, the accreting gas is so tenuous that it cannot cool effi- 
ciently, the viscous energy is stored in the protons as thermal 
energy and eventually advected onto the nuclear compact ob- 
ject. The ADAFs are therefore characterized by small radia- 
tive efficiency and accretion rates (m = M /Msdd < 10 -16 ; 
Narayan & Yi 1995; Rees et al. 1982). At low m the hard X-ray 
emission is mainly due to bremsstrahlung emission from a pop- 
ulation of ^100 keV electrons, and is therefore much harder 
than typically observed in Seyfert galaxies (Nandra et al. 1997; 
Turner et al. 1997; Matt 2000). 

Recently, hard X-ray tails have been discovered in the 
ASCA spectra of several elliptical galaxies (Allen et al 2000). 
Photon indices are in the range 0.6-1.5, therefore remarkably 
consistent with the expectations of the ADAF scenario. Unfor- 
tunately, the ASCA energy bandpass is limited to 9-10 keV, 
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and the detection of these hard tails is therefore difficult and 
somewhat model-dependent. 

One of the most intriguing of these ADAF-candidates is 
the nearby (z = 0.0049) elliptical galaxy NGC 1052. It is 
a narrow-line LINER (Heckman 1980), with a compact ra- 
dio core (diameter ~0.14 pc) and a radio halo of about 3 kpc 
diameter. NGC 1052 is the first type 2 LINER where broad 
lines in spectropolarimetric measurements have been discov- 
ered (Barth et al. 1999). The ASCA spectrum (Guainazzi & 
Antonelli 1999, G99; Weaver et al. 1999) is extremely flat: 
a formal fit with a simple power law yields a photon index 
r ~ 0.1. This suggests either high photoelectric absorption 
(column density /V H — 10 23 cm -2 ) or a Compton-reflection 
dominated spectrum (which, of course, would imply an even 
more strongly absorbed active nucleus). In the former scenario 
the intrinsic spectral index is still very flat (r ~ 1.2). This is 
intriguingly consistent with the ADAF scenario. 

The ASCA results prompted a BeppoSAX program of 
observations of narrow-line LINERs with broad spectropo- 
larimetric lines, mainly focused to investigate their hard X- 
ray spectrum (>10 keV). The results of this experiment on 
NGC 1052 are reported in this paper. 

2. Observation and data reduction 

BeppoSAX observed NGC 1052 between 2000 January 11 
(19:27 UT) and 13 (08:57 UT). The instruments were oper- 
ating in nominal direct modes. Data reduction and analysis 
followed standard procedures, as detailed, e.g., in Guainazzi 
et al. (1999). Scientific products from the imaging instru- 
ments (Low Energy Concentrator Spectrometer, LECS; Par- 
mar et al. 1997; Medium Energy Concentrator Spectrometer, 
MECS; Boella et al. 1997) were extracted from circular re- 
gions around the galaxy X-ray centroid of radii 2' and 4', re- 
spectively. Background subtraction was performed using deep 
blank sky field exposures, accumulated by the BeppoSAX Sci- 
ence Data Center (SDC). The background-subtracted net count 
rates are (1.26±0.09) x 10~ 2 s" 1 and (3.71±0.08) x 10~ 2 s _ \ 
in the LECS (0.1^ keV) and MECS (1.8-10.5 keV), respec- 
tively, corresponding to exposure times of 25.8 and 63.2 ks. 
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Fig. 1. Spectrum plus best-fit model (upper panel) and resid- 
uals in units of standard deviations [lower panel) when the 
model "P" is applied to the broadband BeppoSAX NGC 1052 
spectrum 

Spectral analysis made use of the latest response matrices re- 
leased by the SDC in January 2000. 

In what follows: energies are quoted in the source rest 
frame, and errors are at 90% confidence level for one inter- 
esting parameter (Ax 2 = 2.71), unless otherwise specified. 

3. Results 

3.1. The PDS spectrum (E > 10 keV) 

The 60.0 ks exposure time on NGC 1052 allowed a PDS detec- 
tion, with a total 13-200 keV count rate of 0.22 ± 0.04 s -1 . 
If we take into account the typical systematic uncertainties 
associated with the PDS background subtraction algorithm 
(Guainazzi & Matteuzzi 1997), the detection is at a level higher 
then 5(7 in the 13-90 keV band (see Fig. m. No known hard X- 
ray bright source is present in the (1.3° r PDS field of view. 
The number of expected sources with flux equal or higher than 
NGC 1052 is < 0.06 according to the Cagnoni et al (1998) 
logN-logS relation. A simple power-law fit of the PDS spec- 
trum yields a very flat index: T = 1.0±q |. 

3.2. The broadband BeppoSAX Spectrum ( 0. 1-100 keV) 

The source did not show any significant X-ray variability dur- 
ing the BeppoSAX observation. The reduced \ 2 , when a con- 
stant line is fit to the At = 5760 s light curve are x 2 = 0.76 and 
X 2 V = 0.78 in the 0.1-2 keV (LECS) and 2-10 keV (MECS) en- 
ergy bands, respectively. We will therefore focus in this Section 
on the time-averaged spectra only. 

A photoelectrically-absorbed single component model pro- 
vides an inadequate fit of the broadband (0.1-100 keV) Bep- 
poSAX spectrum {e.g.: x 2 — 185.4/80 degrees of freedom, 
dof, if a power-law model is employed). On the other hand, 
a very good fit (x 2 — 1.02-1.03) is obtained with a two- 
component model, constituted by an absorbed (Ah — a few 



10 23 cm -2 ) power-law plus a "soft excess" below 2 keV. The 
limited statistics prevents us from unambiguously characteriz- 
ing the latter component. In the following we will discuss, as il- 
lustrative examples, models where this soft excess is described: 
either with a power-law ("P" model hereinafter), whose index 
is held fixed to that of the high-energy absorbed power-law 
(thus modeling reflection of the nuclear continuum, scattered 
by an electron plasma - "warm mirror" - along our line of sight; 
Antonucci & Miller 1985); or with thermal emission from a 
collisionally ionized, optically thin plasma (mekal model in 
Xspec; "M" model hereinafter). Model "P" also describes a 
geometry, in which the absorber only covers a fraction of the 
line of sight. The best-fit parameters are reported in the up- 
per panel of Tab. |l[ As already suggested by the analysis of 
the PDS spectrum alone, the absorbed power-law component 
is rather flat (T ~ 1.4). If the power-law in model "P" is sub- 
stituted by a thermal bremsstrahlung (x 2 = 83.6/77 dof), its 
temperature is 150±[qq keV. 

The addition of a narrow Gaussian emission line is required 
at the 98.9% confidence level, according to the F-test, in the 
"P" model (Ax 2 |p = 10 for a decrease of the degrees of 
freedom by two), whereas only at the 90.9% level in the "M" 
model (A% 2 |m = 5). The Gaussian line centroid energy is 
consistent, within the statistical uncertainties, with K Q fluores- 
cent emission from neutral iron. However, the EW of the iron 
line system is too large to be produced in transmission by the 
same cold matter, which is responsible for the attenuation of 
the X-ray continuum (which would imply EW ~ 130 eV for 
a spherical distribution of matter; Leahy & Creighton 1993). 
No iron emission line is expected from an ADAF. The slight 
difference in EW between models "P" and "M" (in the lat- 
ter the iron line profile is partly accounted by the emission 
of the thermal plasma) may suggest a multi-component struc- 
ture of the iron line, which is unresolved by the MECS. Ion- 
ized iron lines could be also produced by the "warm mirror" 
(Netzer & Turner 1997). We have therefore repeated the fit 
in the "P" scenario, assuming that the iron emission actually 
consists of two components: one neutral (E c — 6.4 keV) and 
one He-like (E c = 6.7 keV). The fit is of comparable qual- 
ity (x 2 = 79.8/77 dof), with: EW(6A keV) = 170±}^ eV; 
EW(6.7 keV) = 180±^g eV. The soft excess continuum flux 
at 6 keV is about 1/3 of that of the transmitted component. 
Therefore, the EW of the ionized iron line against its proper 
continuum would be of the correct order of magnitude if pro- 
duced in a "warm mirror" (Matt et al. 1996). The neutral com- 
ponent EW is now consistent with being produced in trans- 
mission by the same matter covering the active nucleus, if its 
covering factor is large. 

A hard X-ray continuum could be in principle due to Comp- 
ton reprocessing of the nuclear continuum, by either the ac- 
cretion disc (George & Fabian 1991; Matt et al. 1992) or the 
molecular torus encompassing the active nucleus (Ghisellini et 
al. 1994; Krolik et al. 1994). This scenario does not, however, 
match our data. A fit, where the absorbed high-energy com- 
ponent is a bare face-on Compton-reflection (model pexriv 
in XSPEC; Magdziarz & Zdziarski 1995) is statistically unac- 
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Table 1. Best-fit parameters and results when the models "M" and "P" (details in text) are applied to the NGC 1052 broadband 
spectrum of BeppoSAX {upper panel), ASCA-ROSAT (after G99; central panel), and ASCA-BeppoSAX-ROSAT {lower panel). 
f s is the scattering fraction (defined as the 2-10 keV flux ratio between the transmitted and the scattered power-law components). 
E c and EW are the centroid energy and the equivalent width of the emission line, respectively. 

"BeppoSAX 
b ASCA 



ceptable (x 2 = 112.6/76 dof). The addition of a Compton- 
reflection component to the "P" model (where only the rela- 
tive normalization between the reflected and the direct com- 
ponent, R, and the intrinsic power-law cut-off energy are left 
free parameters in the fit; an inclination angle of 30° and solar 
abundances are assumed) does not significantly improve the fit 
(X 2 = 79.1/75 dof). The 90% upper limit for two interesting 
parameter on R is 0.6 (for T < 1.65). These results allow us to 
rule out one class of models, which adequately fit the ASCA- 
ROSAT spectra, hence favoring the G99 transmission scenario. 

The observed flux in the 0.5-2 keV (2-10 keV) energy band 
is 0.4 (4.0)xl0~ 12 erg cm~ 2 s _1 . This corresponds to a lumi- 
nosity of 0.4 (4.2) xlO 41 erg s -1 . 

3.3. Comparison with ROSAT and ASCA data 

In the central panel of Tab. |[ the best-fit parameters are re- 
ported, when the "M" and "P" models are applied to the ASCA 
and ROSAT spectra of NGC1052 (see the Tab. 1 in G99). 
NGC 1052 was comparatively bright during the August 1996 
ASCA (2-10 keV flux -3.5 x 10~ n erg cm" 2 s" 1 ) and the 
January 2000 BeppoSAX observations. The only spectral pa- 
rameter showing a significant difference is the absorbing col- 
umn density, which was about a factor of 2 higher in the 
later BeppoSAX observation. The spectral indices measured 
by BeppoSAX tend also to be slightly softer, but still consis- 
tent with the ASCA-ROSAT measurements within the statisti- 
cal uncertainties. 

We have performed a simultaneous fit of the ROSAT, 
ASCA and BeppoSAX spectra, to check whether the improved 
statistics allows us to distinguish between the "M" and "P" 
models. In these fits only the column density absorbing the 
primary nuclear continuum has been allowed to vary indepen- 
dently in the BeppoSAX and ASCA-ROSAT models (ROSAT 
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Fig. 2. Spectral index versus column density iso-x 2 contour 
plots, when the ""P" model is simultaneously applied to the 
ASCA, BeppoSAX and ROSAT spectra. The column den- 
sity is left free to vary independently in the BeppoSAX and 
ASCA/ROSAT models {labels). The contours correspond to the 
68%, 90% {solid lines) and 99% {dashed line) confidence lev- 
els for two interesting parameters. 



spectra are basically insensitive to column densities of the 
order of 10 23 cm -2 ). Normalization constants have been in- 
cluded as free parameters in the models, to account for the 
different fluxes measured in the three observations. The re- 
sults are reported in the lower panel of Tab. [j]. The two models 
yield comparably good fits (x 2 |p = 386.0/445 dof; x 2 |a/ = 
383.0/444 dof). Better data quality is needed to resolve this is- 
sue. The spectral index is indeed much better constrained than 
by BeppoSAX data alone, and still very flat (see Fig. fy. 
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Fig. 3. Spectral index versus extrapolated 2-10 keV flux iso-x 2 
contours if a simple power-law model is applied to the PDS 
spectrum of NGC 1052 and of the "flat" Seyfert 2s of Turner 
et al. (1997). The levels are 68% (solid line) and 90% (dashed 
line) for two interesting parameters 

4. Discussion 

NGC 1052 is the first "type 2" LINER, where broad lines were 
discovered in spectropolarimetric measurement (Barth et al. 
1999). One may therefore expect that its X-ray spectrum re- 
sembles that of Seyfert 2 galaxies, where the luminous nuclear 
emission is photoelectrically absorbed by matter with a sub- 
stantial column density. The observation performed with Bep- 
poSAX demonstrates that this is indeed the case, thanks to the 
first measure of the broadband (0.1-100 keV) X-ray spectrum 
of this object. 

Even more intriguing is the nature of the accretion flow in 
this AGN. The nuclear spectrum is rather flat (T ~ 1.4), as the 
ASCA observation had already suggested (G99; Weaver et al. 
1999). The best-fit value is inconsistent at 2a level with the dis- 
tribution of spectral indices observed in Seyfert Is ((Tsyi) = 
1.87; o-syi = 0.22; Nandra et al. 1997; Reynolds 1997) and 
radio-quiet quasars ((Fr,qq) = 1.93; ctrqq = 0.22; Reeves & 
Turner 2000). It is closer to the distribution of spectral indices 
in radio-loud quasars ((Trlq) = 1-6; ctrlq = 0.15; Reeves 
& Turner 2000). However, NGC 1052 radio-loudness, accord- 
ing to the Wilkes & Elvis (1987) definition, is 1.1, therefore 
pointing to a radio-quiet or, at most, a borderline object. There 
are some Seyfert 2s, whose ASCA spectral indices (Turner 
et al. 1997) are as flat as those measured in NGC 1052. Re- 
cently, Malaguti et al. (1999) pointed out that at least in one 
case (NGC 2110) this could be due to an incorrect modeling 
of a complex photoelectric absorber in the ASCA bandpass. To 
verify this hypothesis, we have extracted from the BeppoSAX 
archive the PDS spectra of the "flat" Seyfert 2s of the Turner et 
al. (1997) sample. In Fig. |we compare the T versus 2-10 keV 
(extrapolated) flux iso-x 2 contours for NGC 1052 and the "flat" 
Seyfert 2s (the results for NGC 5252 are not shown, because 
the PDS detection is too weak). In all cases the intrinsic spec- 
tral index as measured by the PDS is significantly steeper than 
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Fig. 4. Radio-to-X-ray SED of NGC 1052 

in ASCA. This suggests that the nature of the accretion flow 
in the nucleus of NGC 1052 is qualitatively different from that 
normally at work in bright nearby AGN. 

The NGC 1052 high-energy spectrum is can be also 
well approximated by a thermal bremsstrahlung with kT ~ 
150 keV, in agreement with the expectations of the ADAF sce- 
nario. No estimate of the nuclear black hole mass is avail- 
able for NGC 1052. If we use the bulge B magnitude (Ho et 
al. 1997) versus black hole mass relation of Magorrian et al. 
(1998), we obtain M B h ~ 10 8±1 M Q . The bolometric luminos- 
ity of the ADAF component extrapolated from the BeppoSAX 
measurement is ~5 x 10 42 erg s _1 , suggesting an accretion rate 
777 ~ 10" 4±1 .This value is significantly lower than the critical 
77i, below which the onset of the ADAF regime would occur 
(Narayan& Yi 1995). 

Critical diagnostic tools for the existence of an ADAF are 
the radio-to-X-ray Spectral Energy Distribution (SED) and lu- 
minosity ratio (Di Matteo et al. 2000). In Fig. ^ we show the 
NGC 1052 radio/X-ray SED. The radio data points are taken 
from the NED archive. We observe a good qualitative agree- 
ment between the shape of the SED in NGC 1052 and in the 
"candidate-ADAF" elliptical galaxies of the Allen et al. (2000) 
sample. In Tab. ^ we compare the ratios between the powers at 
6 keV and 5 GHz in these galaxies and in NGC 1052. In the lat- 
ter this ratio is well within the (rather narrow) range observed 
in the former objects, confirming the qualitative similarity. For 
comparison, the value of this ratio is typically > 10 5 in radio- 
quiet quasars (Elvis et al. 1994) 

It is also interesting to compare the X-ray with the IR 
photometry of Becklin et al. (1982). Their 1.2/^m flux in 
the innermost 2" (~300 pc) - 26 mJy - corresponds to a 
1/(1.25 fim)/L(5 keV) ratio of about 15. Such a high value 
is typical of Seyfert 2 galaxies (Mass-Hesse et al. 1995). The 
IR emission is therefore likely to be dominated by heated dust, 
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Table 2. X-ray to radio power ratios in NGC 1052 and the can- 
didate ADAF elliptical galaxies of Allen et al. (2000). The data 
for these objects are taken from Di Matteo et al. (2000) 

totally masking the ADAF IR intrinsic emission and therefore 
preventing us from checking the presence of a wind in the 
ADAF scenario (Di Matteo et al. 2000). 

The luminosity of the broad component of H^ is ^2 x 
10 41 erg s" 1 (Ho et al. 1997; Barth et al. 1999). An ADAF can- 
not provide the corresponding amount of (mainly UV) ionizing 
photons. This points to a flow with a small sonic radius r s , con- 
sistent with models where the viscosity parameter a is ^ 10 -2 
(r s ~ a few R g ; Narayan et al. 1997). Outside this radius the 
flow would smoothly connect to a standard thin accretion disk. 

The BeppoSAX data alone cannot rule out the possibility 
that a substantial contribution to the hard X-ray emission comes 
from a nuclear relativistic jet. Chandra will have enough spa- 
tial resolution to image X-ray structures on scales as small as 
~100 pc. However, it is to be noted that the ROSAT HRI im- 
ages of NGC 1052 unveiled only a soft X-ray extended emis- 
sion on a much wider scale (~15 kpc) than that of the radio 
structures (G99), therefore probably due to the emission of dif- 
fuse intergalactic gas. The contribution of this component to 
the hard X-ray emission is most likely negligible. 
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